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Abstract 

We study exact self duality in the model of a Brownian particle in 
a washboard (WB) potential which describes a Josephson Junc- 
tion (JJ) coupled to an environment, for arbitrary temperature 
and arbitrary form of the spectral density of the environment. 
To this end we introduce the Quantum Dissipative Villain Model 
(QDVM), which models tunneling of a degree of freedom cou- 
pled to a linear quantum environment through an infinite set of 
states. We derive general exact mappings on various dual dis- 
crete representations (one-dimensional Coulomb gases or surface 
roughening models) which are exactly self-dual. Then we show 
how the QDVM maps exactly onto the WB model and use dual- 
ity relations to calculate the leading terms of the total impedance 
of a J J circuit, for general frequency dependence of the spectral 
density of the environment and arbitrary temperature. 
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1 INTRODUCTION 



Tunneling in quantum dissipative systems can be described in terms of one 
(or few) special degree of freedom coupled to an incoherent environment of 
quantum oscillators [Q. If the special variable has an underlying discrete 
character the model may have various discrete dual representations, each 
corresponding to specific dressed excitations. 

As an example of underlying discrete character of a continuous degree of 
freedom we may consider a Brownian particle in a double- well potential 
a model which also describes a mesoscopic SQUID 0. In some region of 
the parameter space the relevant physics is well represented by a two-state 
approximation, each state being related to a minimum of the potential. 

In this paper we will focus on the WB model || which describes the 
Brownian motion of a quantum particle in a tilted washboard potential and 
the quantum dynamics of mesoscopic J J circuits 0, For JJs the above 
special variable is the phase if across the junction (see fig.l) a collective de- 
gree which interacts whith all the microscopic variables entering the problem 
(quasiparticles in the electrodes Q, excitations of the circuit 0) modeled by 
an environment with given spectral density. The partition function is 



where m is the mass and a(u) is the damping kernel for the Brownian parti- 
cle [|l|. They correspond to C and Z(uS) in the JJ circuit of fig.l, Zt{oj) being 
the impedance seen by the junction || [7|, ||, Zt(uj) = (1/Z(u) + iuoC) 1 . 
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The WB model accounts for a variety of effects (eg. Josephson effect, 
phase diffusion, macroscopic quantum phenomena, Coulomb blockade, Bloch 
oscillations) in different regimes [|J , due to its "dual" structure. Indeed ip and 
the charge Q at the junction are conjugate quantum variables, [hip/2e, Q] = 
ih and both have an underlying discrete character (it comes from the ex- 
istence of minima in the WB potential for tp, whereas for Q it is related 
to the fact that Q varies in units of 2e due to Cooper pair tunneling). As 
a result pairs of "dual" effects are observed in J J circuits |4j], for instance 
the Josephson effect (a finite current flowing at zero voltage, corresponding 
to classical dynamics of ip) and the Coulomb blockade (finite voltage drop 
with zero current, corresponding to classical stochastic dynamics of Q). It 
is tempting to think that "quantum duality" implies a phase transition trig- 
gered by the ratio V/E c , E c = e 2 /2C being the charging energy ||, but 
this is wrong: the low-frequency response of a single junction device in a 
circuit is determined solely by the properties of the coupling with the envi- 
ronment ||. For instance, for an ohmic environment the critical parameter 
is a = Rq/Z(u = 0) (Rq = h/4e 2 ~ 6.5 kfl is the quantum of resistance): 
Coulomb blockade is possible if a -C 1 (effectively current biased junction), 
whereas for a ^> 1 (effectively voltage biased junction) the dynamics of ip is 
semiclassical. 

In general, collective excitations and duality are more involved than the 
Kramers- Wannier ip-Q duality because they depend mainly on the environ- 
ment || 0. Duality in the WB model with ohmic dissipation at T = 
was studied first by Schmid who found a low-frequency approximate self- 
duality with a 1/a. A dual mapping onto a one-dimensional Tight- 
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Binding model was then introduced [[| for arbitrary temperature, which was 
then used by Sassetti et al. [|10j to map into each other the sub-ohmic and the 
super-ohmic environment. Very recently it was shown that the Schmid self- 
duality is exact at T = for a strictly ohmic environment [|TI|] the m, C — > 
and Zt{oj) = R limit of eq.(||]), where the problem is exactly solvable 

In this paper we study the discrete representations and duality for the 
WB model. We introduce in sec.0 the QDVM and its discrete dual and self- 



dual representations, surface roughening and Coulomb gas models fL3|| . The 
QDVM, besides illustrating in a simple but exact mathematical framework 
how duality works in connection with discreteness and dissipation, can be 
reduced exactly to the WB model with a given environment (sec.[3]). thus, 
also the WB model has an exact self- dual structure for arbitrary environment 
and temperature. As an application we finally calculate the impedance of a 
Josephson circuit for low frequencies. 

2 THE QDVM 
2.1 The model 

The QDVM is defined on the discretized (slice 1/A) imaginary time axis as 



Z(J) 



~[dip T 



. e~5S E T V( VT -2nm. T ) 2 _ ^ 

{m T } 

The discrete variable m T is coupled to <p T which undergoes gaussian fluctua- 
tions tuned by the kernel A® T ,, J T being the source. We can interpret eq.(||) 
as the discretized path integral for <f(r) — 27rm(r), a continuous variable with 
an underlying discrete character (see fig. 2a) coupled to an environment. 



By comparing eq.([3|) with eq.(|T]) it is apparent that the QDVM can be 
obtained from the discretized version of the WB model by using the Vil- 



although it is not valid in the continuum limit V/A — > 0. Notice however 
that here we use a strategy different from ref. we first discuss the dual 
properties of the QDVM and then we show how to obtain exactly from it the 
WB model, in the limit V/A -> 0. 

2.2 Discrete representations of the QDVM 

Surface roughening representation Starting from eq.(|3|), we can elimi- 
nate {(fi T } and obtain a surface roughening model [13|] on the imaginary time 
axis, where m T is the surface height in [r, r + 1/A] (see fig. 2b). Each element 
m T of the surface interacts with all the other elements. 

Coulomb gas representation I: e charges We now introduce the charges 
e T = m T+ i/^—m T and perform a discrete double integration of the interaction 
to obtain the generating functional of the QDVM in the e-representation 



This is a set of interacting (A TT /) integer charges e T g] — oo, oo[ (see fig. 2c). 
Z v {J) is the generating functional for the damped harmonic oscillator 
and the interaction is given by 



lain approximation |TJ[], so we can identify the kernel A® T r. The Villain 
approximation has been successfully used to study mesoscopic J J arrays |T5[ , 



K} 



e T A TT , e r ,+ £ T J^e T 



(4) 




{4TT 2 m-47r 2 a{u)/2u 2 } V 
muj 2 - a(uj)/2 + V 
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The source jf^ can be obtained from the source J T of the QDVM \W . 

If we classify the configurations in eq.(f|) according to the number of 
charges of each species we obtain the Coulomb gas representation 



2 (e) (J (e) ) = E 

{Me} 



n wr 1 



{n} 



Here each term of the summation is labeled by a set {.Me} where Ai s is the 
number of charges e^O present. The charges are numbered with the index 
/ = 1, . . . , Mt, where Mr = X^-oo M s - For each set {M s } only a special 
configuration enters eq.(||). The properties of the kernel enforce the charge 
neutrality condition J2gL-oo -Mg-e = 0. 

Coulomb gas representation II: n charges A different Coulomb gas 
representation can be obtained if we apply the Poisson transformation to 
eq.(||) (or to eq.(||)). We obtain the n-representation of the QDVM 

Z{J) = Z°(J)-ZW(JW) 
Z^ n \j^) = V ^y= n - v tt> n T'+T, T Jr n) n T ^ ^ 

{nr} 

which is again a gas of interacting integer charges n T e] — oo,oo[. The 
partition function Z Q (J) of the free Brownian particle ]]] appears and again 
can be obtained explicitly from J T fl6|| . The generating functional 
Z^ n \j'^) has the same form of eq.(|]) with a different interaction given by 

where is the discrete Green's function of the free Brownian particle JT]]. 
By proceeding as before, a Coulomb gas form for Z^^J^) is obtained, 
structurally identical to eq.(||). 
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2.3 Self-duality 

The fact that Z^ e \j^) and Z^ n, \j'^) are structurally identical indicates 
that the QDVM possesses a self-dual structure. The relation between the 
interactions A T and T> T assumes a simple form if we introduce the func- 
tion h°(w) := Z T (-iw)/R Q . Then (1/A) V w = 27i7i (|w|)/|w| + 1/V =: 
(2tv/\uj\) h(\cu\) and (1/A) A w = (2tt/\u\) [l/h(\u\)]. So the exact self dual- 
ity in the QDVM is given by h(u) «-> l/h(u), which is a sort of a Norton 
transformation. 

For low frequencies we obtain the a 1/a and the sub-ohmic <-» super- 
ohmic dualities analogous to the ones possessed by the WB model H, [Ufl. 



Self-duality is a powerful tool for obtaining nonperturbative results. For 
instance from our low-energy results we can infer the existence of the criti- 
cal line a = 1 for ohmic dissipation in the QDVM, analogously to the WB 
model |§ . For any dual mapping, perturbative results in a certain representa- 
tion may correspond to nonperturbative results in a different representation 
(see the next section). An appealing feature of self-duality is that one can 
write exact relations (as the eqs. (§) and (|10D we will use below) and exact 
equations for the correlation functions [pX 



3 THE WB MODEL 

3.1 Mapping of the QDVM onto the WB model 

Now we show that if we take the continuum limit in a special way the 
QDVM on a lattice reduces exactly to the WB model. Indeed from the n- 
representation we obtain the exact perturbation series in V as given in ref . |§ . 
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The proof goes as follows. First one considers the Coulomb gas representation 
for n-charges of the QDVM with the substitution V — > V = A/21n(2A/V). 
Then one shows that the configurations where only the charges n\ = ±1 ap- 
pear give the exact perturbation expansion in V for the WB potential M. 
Finally one shows that all the other charge configurations carry an addi- 
tional factor A with a negative power and give a vanishing contribution to 
the partition function as A — > oo. 

The consequence of this mapping is that the WB model possesses an exact 
self duality for arbitrary T and arbitrary spectral densities which generalizes 
the approximate Schmid self-duality 0. 

Notice that the exact self-dual mapping is not evident if one starts directly 
from the WB model. In fact the perturbation expansion is recasted in a 
Coulomb gas with only ni = ±1 charges, whereas in the instanton expansion 
(the dual Coulomb gas) all kind of charges ti g] — oo, oof may appear. Self- 
duality can be exploited by carrying over the exact results for the correlation 
functions of the QDVM to the WB model. 

3.2 Impedance of a J J circuit 

As an application we calculate the leading terms of the retarded (tptp) corre- 
lator which is related to the impedance Z(f2) seen by the source of the circuit 
in fig.l Starting from the thermal ((p(p) correlator of the Villain model we 
perform an analytic continuation, uj n — > —ioo), and then take the contin- 
uum limit with the prescription of sec. |3.1| . It would be desirable to perform 
the analysis in real-time ffl], but at present we content ourself to check the 



QDVM method against available exact results [12]. We we start with the 



8 



exact expression for the correlator (1/A) ((p(p) UJn vs. (nn} Un 



AGlf A(nn) t 



2 r 



4H- 4M " 4">(o)-x 



.£>< 



(8) 



(n) 



where is the retarded correlator of the free Brownian particle, and 
can be evaluated by analytic continuation of the first term of the Coulomb 
gas expansion 



Jo 



Here Fq(s) accounts for free circuit fluctuations [^|, |7|, [L7|, |T| 



LJ R c 



ctgh — (coso>s 



r 



i sm cjs 



In particular QXj^ (u) = (V/2) [-Fo(^) — -Po( — ^)] ; and its Hilbert trans- 
form gives JJXr. We introduced Pq{oj) = J^ds exp{zu;s} exp{F > (s)} 
which can be calculated numerically for almost arbitrary environment at all 

As a check of the theory we calculated analytically the low-frequency 
behavior of Z(Q) at T = 0, using the toy ohmic environment $IZt(uj) / Rq = 
a^ 1 exp{— | uj | /a)}. For a < 1 we find 



9ftZ(fl) 
Rq 



{l-AlT 2+2 / a -M] 2 } (9) 



which reproduces known exact results [fCql , i.e. the leading behavior of the 
correlation function is ~ f2 2 for < a < 1/2 and ~ f2~ 2+2 / a for 1/2 < 
a < 1. It corresponds to a vanishing junction admittance, which describes 
the Coulomb blockade. The coefficients in (||) depend on the details of the 
environment at frequencies of the order of the cutoff uo and can be evaluated 
numerically. 
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Then we start with the exact expression of the (<p<p) correlation function 
in the e-representation and then perform the analytic continuation 

2 




UJ 



X%{Q)-X${ U ) (10) 



where is the discrete Green's function and x\ is the retarded correlation 
function for the damped oscillator. As before we consider the first contribu- 
tion in the Coulomb gas expansion eq.(||) for A(ee) UJn 

2A y e A„-A„ r „ ff -n^zr _ 4A 2 [°° ds e™ Se F> W (11) 

A different correlator for the fluctuations appears, F > (s), coming from the 
analytic continuation of A r , from which we can define a P(u) function which 
enters the expression for Z(f2). To evaluate Z(f2) at all frequencies we need 
a regularization scheme, which is not needed for checking analytically the 
power-low behavior at low frequencies. For a > 1 we find 

Rq 

which in this case gives also the junction impedance. Again known exact 



results JT2| are reproduced. 



4 CONCLUSIONS 

In summary we have exploited the Schmid |§ self-duality in the WB model. 
Using the QDVM we have shown that it is exact, that it holds for finite 
temperatures and that it can be extended to an environment with general 
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frequency dependence of the spectral density. Exact low-frequency results 
for the strictly ohmic environment are reproduced using this method. 

The study of a general environment is important for design and perfor- 
mances of the devices, for instance for an optimal control of the noise. In the 
Coulomb blockade regime it is possible to treat a general electromagnetic 
environment in realistic situations using the P-function theory |5], [5], |17| . 
The QDVM theory can be combined with the P-function theory to calculate 
quantities of direct experimental interest different regimes. 
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Figure captions 



Figure 1: The JJ circuit described by the WB model: the special variable 
is the phase across the junction ip and the impedance Z{uS) is modeled by an 
environment of quantum oscillators. Classical and quantum fluctuations of 
tp are modulated by the fluctuations of the environment. 

Figure 2: (a) A typical trajectory in imaginary time for the coordinate of 
the QDVM, y?(V) — 27rm(r). (b) A typical configuration after the elimination 
of {tp T }'- the set {m T } is associated to a one-dimensional interface; each 
elementary (1/A) piece of surface interacts with all the other elements, (c) 
The corresponding charge configuration in the Coulomb gas representation 
eq.(||): in general charges may take any integer value in e/ g] — oo, oof. 
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Falci and Weiss - Figure 1 
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Falci and Weiss - Figure 2 
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